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Abstract: We present a phenomenological appraisal of the results of several searches for Supersym- 
metry (SUSY) performed at the LHC by the CMS collaboration and interpreted in the context of 
the cMSSM. Part of the analysis focuses on which SUSY production processes are being probed. We 
observe that much of the current exclusion region is dominated by squark-squark and squark-gluino 
production, and explain the shape of the exclusion curves. Based on this analysis and an estima- 
tion of future reach, additional simplified models are suggested. Other phenomenological details are 
discussed, such as the effect of radiation on acceptance. 
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1 Introduction 

Significant effort at the LHC is dedicated to anaiyzing data to find evidence of, or to exclude, models 
of Supersymmetry (SUSY) (Refs. [1, 2] provide thorough theoretical reviews.) Many of the results 
of these analyses are interpreted in the context of the constrained Minimal Supersymmetric Standard 
Model (cMSSM) [3], with exclusion contours shown in the plane of the universal scalar mass toq and 
the gaugino mass mi parameters, with other cMSSM parameters fixed (tan/3, s(7n(/z)). A recent 
result from CMS [4] is presented in Fig. 1, based on 1.1 fb"-'^ of data. 

Our initial attention is drawn to the curves in Fig. 1 labeled Jets+mHT [5], ax [6], and MT2 [7], 
referring to limits from specific analyses known by those colloquial names. These analyses are sensitive 
mainly to the production and hadronic decays of sparticles, whereas the other analyses, with limit 
contours labeled 1 Lepton, SS/ OS Dilepton\^-lQ] , depend explicitly on branching fractions of sparticles 
to leptons (either directly or through decays of gauge bosons) . 

The hadronic limits demonstrate that: 

• for low mo, there is a significant reach in mi (~ 540 GeV); 

• for large moi mi is limited to smaller values (^ 300 GeV). 

The parameters mg and mi are not the physical masses of the theory, which set the size of production 
cross sections and the kinematics of decays. However, is basically set by mi {Mg ^ 2.5m i), while 
Mq^ has a quadratic dependence on mo and mi . From Fig. 1, it is clear that the excluded gluino mass 
is not a flat line: the exclusion depends strongly on the squark mass. (Figure 18 of the Appendix shows 
the correlations between Mg and Mg in the (mo, mi) plane relevant for this analysis.) In particular, 
one can see from the contours of fixed gluino and squark mass, labeled q and g, that Mg < 1 TeV is 
excluded regardless of Mg, whereas Mg < 1 TeV is only excluded for Mq <1.1 TeV. The exclusion 
curve seems to follow the semi-circular contours of fixed Mg for low mo < 1 TeV, and the flat contours 
of flxed Mg for larger mo. 
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Figure 1. Observed limits from several 2011 CMS SUSY searches plotted in the cMSSM (mo,mi) plane, 
with tan /? = 10, Aq = 0, sgn(p) > 0. 



The physical masses of the cMSSM are highly correlated, as expected from a model that has 
many observables but only a few parameters. While the quantitative structures of the cMSSM depend 
on the exact specification of mQ,mi, Aq, tan l3, sgn(ii), the qualitative features are set by mo and 
nil (for larger values of tan/3, third generation sparticles may be lighter. From the collider physics 
phenomenology point of view, this affects the other spectra mildly, but can have significant implications 
for decay modes.) It is well known that Mq > Mg in most of this plane and that the gaugino masses 
come in specific ratios (M^^^ = M^o ^ One can question whether general results can be 

learned from examining the limits within the cMSSM, and for a single choice of cMSSM parameters. 
In particular, this limit is for tan/3 — 10, Ao = GeV,sgn{iJ,) > 0. (We also comment on the case 
of tan/3 = 40,^0 = —500 GeV, sgn{fi) > 0.) We argue that several lessons can be learned from this 
particular cMSSM interpretation, despite the fact that it is explicitly model-dependent. 

In this paper, we will investigate: 

1. the dominant SUSY processes that contribute to these limits. In particular, we would like to 
understand the shape of the CMS exclusion curves in Fig. 1. Analyzing the plane this way is a 
mapping of a point in (mo,TOi) into simple or simplified models; 

2. the impact of this mapping on the selection of simplified models for future studies to better 
understand the experimental results; 

3. the sensitivity of the various analyses to specific SUSY production processes and decay modes; 
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4. phenomcnological details relevant to this analysis and other new particle searches, such as kine- 
matic features of various processes and the role of radiation in searches. 



2 Dominant Processes 

The purpose of this section is to examine which processes are contributing to our current exclusion 
contours in the cMSSM. Here, we concentrate on the size of production cross sections without con- 
sideration of a specific selection of events, as in a search analysis. We will justify this later. The 
deconstruction of the cMSSM into subprocesses is a mapping into a set of simplified model spectra 
(SMS) (see [17], for example). This will help us to understand the generality of these cMSSM results, 
and interpret them in simplified models. The main result is presented in Fig. 2, which shows the 
fractional importance of various SUSY pair production processes as a function of (mo, mi) for the 
same parameters as in Fig. 1. 
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Figure 2. Fractional size of various subprocesses as a function of (mo, mi) for tan/3 = 10, Aq = 
GeV, sgn{fi) > 0. 



The labels in Fig. 2 represent (from top-left to bottom-right) the following production processes: 
sg = squark-gluino, ss — squark-squark, tb = stop-antistop, nn — neutralino/chargino pair, ns = 
neutralino/chargino-squark, sb ~ squark-antisquark, gg — gluino pair, II — slepton-antislepton, ng = 
neutralino/chargino-gluino. The classifications sum over similar processes: nn, for example, includes 
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all NiNj, CiCj, NiCj processes. (In practice, only a few processes tend to saturate the sum, e.g. C1N2 
and CiCi in this case.) For brevity, the mostly-negligible sbottom-antisbottom pair production is not 
shown in Fig. 2. These and other production cross sections were calculated using the computer code 
Prospino to NLO accuracy [11]. The physical particle properties were calculated using the MSSM 
evolution code SoftSUSY [12], while sparticle and Higgs boson decays widths were calculated from 
SUSYHIT [13] (sparticle decays are considered later). 

The all-hadronic searches naively target squark and gluino production. Figure 2 demonstrates 
that squark-squark production and squark-gluino production have large rates at small mp. For small 
mi, gluino pair production is important, while the production of electroweak gauginos is important, 
if not dominant, over most of the plane. 

For further clarity, we focus on the largest process at each point in the plane, which is presented 
in Fig. 3. Only four classes of processes every satisfy the criterion of being the largest cross section 
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Figure 3. The largest cross section in the (mo, mi) plane, with contours showing the fractional size. 

at any point (a small region of slepton production for small toq and large mi is not shown). When 
squark-squark production is the largest (roughly 400 < mo < 500 GeV, 100 < mi < 300 GeV), it 
comprises 30-40% of the total SuSY rate at that point. ^ Gluino-pair production is largest in only 

^The large contribution of ss processes to the total SUSY production cross section was noted for the LMl benchmark 
point in Ref. [20]. 
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a limited region of mi (roughly mi < 120 GeV); the other gauginos have much more significant 
production rates. Clearly, searches that are sensitive to the electroweak gaugino processes (mainly 
C^Ci^ , C^N2) can gain inroads to the upper-right portion of the plane, sg is competitive with ss, and 
the largest process for a the region mo < 400 GeV, 100 < mi < 300 GeV. We reproduce Fig. 2 for 
the case of large tan/3, namely tan/3 = 40, = —500 GeV , sgn{fi) > 0, in Fig. 20 of the Appendix. 
The qualitative behavior is the same. 

These conclusions are noteworthy. First, squark-squark production can only occur if a gluino 
exists, whereas squark-antisquark production arises directly from the QCD Lagrangian. Even if the 
gluino is heavy enough to have a suppressed production cross section, the squark-squark production 
cross section can be substantial (we study this more quantitatively later). Thus, an indirect way 
to prove or disprove the existence of a gluino with a TeV-scale mass is to establish that a signal is 
squark-squark or squark-antisquark production. Potential clues could be obtained through jet-charge 
or exploiting some differences in gluon radiation arising from the contributions of valence or sea quarks 
to the process (see later discussion). 

Squark-gluino production, squark-squark, and gluino-gluino production processes are all compet- 
ing at once. We will investigate what mass relations between the squark and gluino allow squark-gluino 
production to dominate over the other possible production processes. 

The nn processes are mainly the SUSY duals of WW and WZ production, which is typical of 
the cMSSM. Almost all of nn is Cf N2 and Cf C^ production, with the lightest chargino Ci and the 
second heaviest neutralino N2 nearly degenerate in mass. 

The cross sections for stop and sbottom pair production are never competitive with other processes 
(even with tan/3 = 40) for these cMSSM parameters. We will discuss this more later. 

The total SUSY production cross section in the (mo, mi) plane is presented in Fig. 4. The 
contours of constant cross section are nearly horizontal lines in mi . This is not surprising, since the 
cross section for the nn process is set mainly by the chargino and neutralino masses. Comparing with 
the limits in Fig. 1 (and reproduced in this figure), we observe that current searches with 1.1 fb~^ of 
data are excluding the cMSSM points that yield a total cross section of several hundred fb, or several 
hundred SuSY events before event selection. 

3 Simplified Models 

The simplified model approach originated as a topological classification of signals to allow hypothesis 
testing [15, 16]. It is based on the observation that the kinematics of new particle production and 
decay is often characterized by the mass scales involved, and not by the particulars of matrix elements. 
In the absence of signals, simplified models (SMS) can be used to present results with less model 
dependence [17]. Thus, SUSY search limits can be recycled to provide information to the community 
about the dependence of these limits on the masses and decay chains in a particular model. SMS are 
also useful to the analysts, since it reveals the coverage of the cuts used in a particular analysis. 

By breaking down the results in the cMSSM into individual processes, we are mapping into 
simplified models. Do we have the right ones or enough? 

This simple analysis of only the production cross sections can already have an impact on the choice 
of simplified models used to interpret the data. The main suggested topologies, Ijg and qq* [19], are 
not the most relevant in most of the cMSSM plane, even though we have relied on these processes to 
motivate our initial choices of topologies. If our aim is to provide motivated SMS that cover the main 
possibilities in a model like the cMSSM, and/or to allow for the deconstruction of a particular point 
in the cMSSM plane in terms of SMS, then other topologies need to be considered, such as: 
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Figure 4. The total SUSY production cross section as a function of (mo, mi) in the cMSSM for tan/3 = 
10, j4o = GeV, sgn{jj,) > 0. The strongest limit curve from Fig. 1 is overlaid with the label Jets+mHT. 

1. squark-squark (ss) production; 

2. squark-gluino (sg) production; 

3. Weak gaugino (nn) production, like Cj^Cf and C^N2- 

However, there are caveats. A specific point in the (mo,TOi) plane will likely not have the simple 
structure of the simplest of SMS, with a single decay chain. Instead, cascade decays like qj^ — ^ 
qN2,qL ^ 1^2, qCi are competing with single step decays like g — > qNi. In general, cascade decays 
tend to soften the spectrum, but can provide more jets that increase acceptance. Thus, limits in 
the SMS may be different than those in the cMSSM for similar choices of Mq, Mg, . 

First, we examine the qq process and compare the kinematics of qq and qq* production. We 
expect some differences from the effect of the parton distribution functions, since valence-valence 
quark scatterings can only contribute to qq production. Differences arise at the hard-process level, 
but also from initial state radiation. 

Figure 3 shows comparisons of the q (or q* ) px after the hard event and parton shower generation 
from Pythia (see the Appendix for details). Comparisons of the rapidity are also shown in the 
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Figure 5. Comparison of the transverse momentum pt of the produced q, q* 
in the Appendix. 



Plots of rapidity y are shown 



Appendix. Values of 150 < Mq < 650 GeV are displayed, though we have studied larger values. In 
general, differences are less pronounced for larger Mg. For Adg < 800 GeV, the pr of the squark is 
harder for ss production, while the rapidity y is less central. Thus, the acceptance of ss events should 
be different from sb. If the acceptance is also dependent upon a jet arising from initial state radiation, 
there will be further differences. 

Since interpretations of data with simplified models have been presented only for gg and qq* 
topologies so far (see, for example, Ref. [21]), we encourage a study of topologies. However, as 
a first estimate, one can ignore the differences in acceptance, and interpret the data using ss cross 
sections, instead of sb. This requires a choice of Mg, since the ss cross section vanishes otherwise. We 
have checked that the kinematic dependence on Mg is negligible, and the acceptance for Mg — 2 TeV 
is nearly identical to that for Mg = 4 TeV, so that signal events can be studied "independently" of 
Mg. Later, we will see that the ss cross section exceeds the sb one even for cases when Mg is quite 
large. 

Next, we turn to the squark-gluino process. We consider the complete MSSM description of 
squarks and gluinos, so that all production cross sections are competing and interconnected. We have 
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discussed already the gg, qq, and qq* contributions. In Figure 6, we show the relative size of each cross 
section as a function of the mass splitting between the gluino and squark, for several fixed values of 
Mq. sg is the dominant production process in a small window, roughly near Mg ~ Mq, that changes 
moderately with Mq. For low Mq, the dominance occurs for nearly degenerate squarks and gluinos. 
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Figure 6. The relative SUSY QCD production cross sections as a function of the mass splitting between the 
gluino and squark. The vertical line is for degenerate Mq,Mg. 



However, for higher Mq, the maximum occurs when the gluino is somewhat lighter than the squarks. 
This leads to a parton-level signature of g(— >■ qqNi)g{-^ QQ^i). Of course, specific cuts might 
select sg over other topologies, so it could be important even what it is not the largest cross section. 
This suggests SMS with Mg within ±200 GeV of Mq. We note also the rise of qq over qq* production 
once Mq > 500 - 600 GeV. For Mq = 600 GeV and Mg = 1200 GeV, the qq production cross section 
is roughly 0.5 pb. 

In Figure 7, we present the same information as a contour plot in the Mq, Mg plane. Note the 
unequal limits for each axes, chosen to emphasize the size of the qq^ cross section for large Mg. Focusing 
on the behavior of the qq cross section at fixed Mq — 400 GeV as a function of Mg, one observes a 
decrease of only one order of magnitude in the cross section as the Mg varies from 275 to 2000 GeV. 

Next, we comment on electroweak gaugino pair production. In the cMSSM, the mass splitting 



- 8 - 



Total Production Cross Section (pb) 



400 £00 SOO 1000 




400 eoo 300 1000 

niQ (GeV) 

Figure 7. Levels of SUSY QCD production cross sections as a function of the gluino and squark masses. 

between the chargino 0.8mi) and the LSP 0.4m i) is roughly 0.4toi . Once mi > ^ 200 

GeV, the decay Ci — >■ WNi is allowed for on-shell particles. From Fig. 1, we observe that the current 
leptonic analyses are already limiting models of this type, since their contours extend up to and beyond 
mi ^ 200 GeV. (We have not studied in detail the shapes of the leptonic limits, but leave this for 
a future study.) This may be viewed as a worst-case in terms of leptonic branching ratio (at specific 
points, N2 — )■ Nih{—>- bb) can occur. We comment on this later). Even hadronic searches may be 
sensitive to these processes, which naively yield 4 hard jets and^T- For heavy enough Ci and N2, the 
gauge bosons produced in decays may be significantly boosted. 

In Figure 8, we show contours of significant decays modes of Ci and N2 for tan /3 — 10, Aq = 
0,sgn{n) > 0. The contour of BR{Ci WNi) ~ 1 is the sofid, black fine near mi ^ 200GeV - 

our naive expectation of where this decay should turn on. The region where N2 — > ZNi is significant 
is small and enclosed by the blue, long-dashed line, and its lower bound is near the start of the 
Ci contour. Near mi ~ 300GeV, a new decay mode is dominant, namely N2 — hNi. This has a 
(negative) impact on multi-lepton searches, but also opens the possibility of a Wkfr signal. The 
dominance of N2 — J> hNi is a feature, or peculiarity, of the correlations inherit to a model like the 
cMSSM, with only a few input parameters. A reasonable extension of the cMSSM is to decouple the 
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Figure 8. Regions wliere BR(Ci -s> WNi), BR(iV2 ZNi), and BR(Af2 hNi) are significant in the 
(mo, mi) plane for tan/3 = 10, Ao — GeV , sgn{n) > 0. 



Higgs scalar from the sfermion scalar parameter. A simple treatment, suitable for collider physics 
phenomenology, would be to remove the Higgs boson from the spectrum entirely, and recalculate the 
decay widths of all sparticles. A consistent approach to this is handled in the PMSSM [18]. 

Based on this, we suggest that a compelling simplified model is CfC^ and C^N2 production, 
with Ci W''*^Ni and N2 ~> Z''*^Ni or N2 ^ hNi. The three-body decays to off-sheU W or Z 
bosons come into play when the mass splitting between mother and daughter is small. = Mj^ 

is a motivated choice, but M^^ should be considered arbitrary. In the CMSSM, there is model 
dependence arising from the q contributions in the matrix-element calculation of the production cross 
section. The q Feynman diagrams are much less significant (numerically) than the Standard Model 
gauge boson ones, which arise naturally from the neutralino/chargino kinetic energy piece in the 
Lagrangian. Reference cross sections can be chosen accordingly, ignoring q contributions. Because 
of the importance of gauge boson diagrams, CfN2 production is roughly twice as large as Cfcf' 
production, and N2 production is roughly 60% larger than N2 production. This mirrors the 
behavior of W* and 7* /Z* in the Standard Model. 

Even though we have downplayed the pair production of heavy-flavor squarks, this view may be 
too model dependent. Contours of stop and sbottom mass for the cMSSM points considered here 
are shown in Fig. 19 of the Appendix. A very large splitting between the stop and/or sbottom and 
the other squark masses is not obtained (except for small mg and rai) for these choices of cMSSM 
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parameters. The stop is special in the MSSM, whether based on fine-tuning arguments or its relation to 
the Higgs boson mass, and one mass eigenstate may be significantly lighter than most other sparticles. 
Later, we observe that the acceptance for tb events is high for generic hadronic selections. Targeted 
analyses, which may leverage the presence of 6 or t quarks in the final state, could have a much larger 
acceptance. Here, we will consider a scenario where Ni , t and g are the only relevant sparticles in 
our current energy regime. Figure 9 shows the relation between gluino and stop pair production as a 
function of mass, where the gluino mass has been offset by the top mass (175 GeV) to emphasize the 
kinematic limit for g — >■ tt. For a stop mass of 300(600) GeV, a gluino of mass 550(950) GeV has a 
comparable production cross section; for smaller mass splittings, gluino production will provide more 
stop events, with additional t quarks in the final state (half with the same sign). 

tb gg 
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Figure 9. A comparison of stop-antistop and gluino pair production as a function of mother mass. The 
gluino mass Mg is offset by rrit = 175 GeV, the kinematic limit for the decay g -> ti. 



4 Studies of the Jets+mHT and ax Analyses 

In this section, we justify our previous statement that an analysis of the production cross section is 
meaningful. To this end, we have performed a generator-level analysis that replicates the Jets+mHT 




and aT event selections. These two analyses are the most straight-forward to reproduce (we leave 
TOT2, which requires a hemisphere decomposition of events and an iteration over invisible particle 
candidates, for a further study). Details of the experimental analyses are to be found in the earlier 
references. There are several results we wish to convey. The first two demonstrate that the current 
analyses are indeed sensitive to the dominant processes, though in different ways. A third shows the 
overlap between the two analyses. The last is to comment on the impact of a leptonic veto on these 
selections. A discussion of the role of extra radiation of quarks and gluons in the analysis results is 
presented in the next section. 

Figure 10 shows the acceptance of events in the {mQ,mi) plane, broken down by process type for 
the Jets+mHT analysis.^ Only six of the process classes are significant. For this plot, the acceptance 
is defined as the number of events that pass the event selection per category, divided by the total 
number of events studies for that value of (mo, mi). The plot demonstrates that the sg, ss, and nn 

Jets + MHT cMSSM Acceptance 
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Figure 10. The acceptance (number of events that pass cuts divided by the total number of events generated 
at the model point) broken down by subprocess type for the Jets+mHT analysis. 

process classes are the most significant after applying the selection cuts. The contours of acceptance 
reflect the contours of total cross section shown in Fig. 2. This implies an acceptance that is relatively 

^We apply the Baseline cuts of [Ht,$t) >{350,200) GeV. See the Appendix for other details. 
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independent of process type, provided that the primary sparticle masses are similar. The nn process 
exhibits a turn-on associated with the threshold in the selection. 

A related quantity is the efficiency of an individual process in isolation from the other processes in 
a full model, like the cMSSM. This efficiency is the relevant one needed for a simplified model analysis. 
Figure 11 shows this unbiased acceptance in the (mo, mi) plane by processes type for the Jets+mHT 
selection. This analysis has a high efficiency for most types of processes, except for the nn class, and gg 
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Figure 11. The acceptance (number of events that pass cuts divided by the total number of events generated 
at the model point for the given process) broken down by subprocess type for the Jets+mHT analysis. The 
contours are drawn in a logj^g scale. 



for lower Mg. There are still correlations built into these results; for example, the ss process depends 
upon the gluino mass. However, the effect on kinematics is less than that on production rates, as 
observed earlier. 

The behavior of the nn, gg, ng acceptance is relatively easy to interpret. All the gaugino masses 
are roughly set by mi independent of mg, and the mass scales sets the hardness of the kinematics, 
and, hence, the selection efficiency. Too low values of , will yield too low to pass these 
selections. The sg acceptance is less obvious. There is some loss of efficiency near the diagonal line. 
From Fig. 18 in the Appendix, we observe that this coincides with Mg ^ Mq. We expect that squark 
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decays to a soft jet and the gluino, and the gluino efficiency is lower because of 3-body decays. 

We have performed a similar analysis for the ar event selection, displayed in Figs. 12 and i3. The 

O.J cMSSM Acceptance 
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Figure 12. The acceptance (number of events that pass cuts divided by the total number of events generated 
at the model point for the given process) broken down by subprocess type for the oit analysis. 

acceptance in the cMSSM is similar, but different, from that in the Jets -/-m^ffT selection. The similarity 
is somewhat surprising, given that ax focuses on the two leading jets, and their relative angle, which is 
characteristic of qq* production.'^ The acceptance per process is shown in Fig. 13. There are patterns 
which are similar to those in Fig. 11, particularly the loss of efficiency near Mg ~ Mq. 

Given the similarity of the exclusion limits (Fig. 1) and the acceptance for the two selections, 
one wonders whether the two analyses are highly correlated, despite the different cuts. In Fig. 14, 
we address the overlap between the two selections, by calculating the acceptance of events that pass 
both selection requirements. Here, we show the fraction of events that pass both selections relative to 
those that pass the Jets+mHT selection. The overlap ranges from 20-50% depending on the process 
and point in the (mo,TOi) plane. There is less overlap for ]fg production and decay, and more for qq 
topologies. This suggests that a combination of these selections would be beneficial. 

handles the case of multiple jets by combining them into two pseudo-jets, whose kinematics are used in the naive 
two-jet formula. This recombination helps recover some jet energy lost to final state radiation. 
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Figure 13. The acceptance (number of events that pass cuts divided by the total number of events generated 
at the model point for the given process) broken down by subprocess type for the ar analysis. The contours 
are drawn in a logjQ scale. 



There is one last issue regarding the hadronic analyses that we wish to address: the impact of 
the lepton veto. A veto is applied on events that have an electron( muon) in the central rapidity 
region of the detector {\ri\ < 2.5) with pT > 10(15) GeV to reduce backgrounds. In Figure 15, we 
show the fraction of events that fail the Jets+mHT selection when the lepton veto is applied at the 
end. The impact of the veto can be severe. The increase near itiq — 900 GeV and mi — 400 GeV 
can be attributed to the turn-on of the decay 5 — )■ t + i, where the t can also decay to a t. As 
we gain confidence in our modeling of Standard Model backgrounds, we may be able to apply more 
sophisticated cuts that reduce the backgrounds without sacrificing signal. This particular example 
also illustrates that simplified models that are too simple may not be able to catch all the complexities 
of a particular model. 



5 Other Phenomenological Considerations 

One may well wonder why production and decay passes an event selection requiring 3 or more jets, 
or why gg production is not constrained at the same level. The presence of additional jet activity can 
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Figure 14. The fraction of events that pass both selections relative to the events that pass the Jets+mHT. 



arise from the evolution of partonic structure through initial or final state radiation (ISR or FSR). The 
probability that a quark produced in the decay of a heavy squark at the scale Qf ^ Mq will remain 
a single parton down to the jet resolution scale is given by the Sudakov form factor, approximately: 

Choosing Mq = 500 GeV and pt ~ 50 GeV, we estimate this probability at roughly 40%, or 16% 
when considering a pair of squarks. Clearly, vetoing a third jet would significantly reduce the signal 
acceptance. The modeling of FSR, while not free from theoretical ambiguity, has been studied in some 
detail at LEP (though for a color singlet system) and in top decays at the Tevatron. The estimate of 
ISR is less under control, because of the dependence on parton distribution functions and the choice 
of phase space for parton emission. Jets from ISR will tend to be uncorrelated with the produced 
sparticles, whereas a selection that includes jets from FSR might recover some of the resonant structure 
in the signal process. 

Figure 16 shows the fraction of events at each (mo, mi) point that pass the Jets+mHT selection 
with an additional jet from ISR. The origin of a jet is traced using Monte Carlo truth information 
and the highest-pT quark or gluon constituent of a jet. Typically 20-30% of the jets arise from ISR. 
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Jets+mHT Lepton Veto 









BSi 




600 




5S0 




50O 












450 








400 








350 


E 






300 




250 




2C0 




ISO 




100 



















1/ 






\ 








/ / / 


































1 1 

2C0 400 


1 1 1 1 1 

GOO ffiO ICCO 1200 1400 

mo (GeV) 


1 1 
leoo isoo 



o.e 



0.5 



■ 0.4 



■ 0.3 



- 0.2 



- 0.1 



0.0 



Figure 15. The fraction of events that fail the Jets+mHT selection because of a veto on identified leptons. 



The region where the gluino is excluded relies more significantly on ISR, presumably because of the 
softening of the jet spectra and-^j- through the 3-body decay g — qqNi. 

The cMSSM has typically large intrinsic ^t, because of the large mass splitting between the q/g 
and A^i. Also, the existence of cascade decays leads to additional hard jets in the final state. This 
ameliorates some of the dependence on ISR. On the other hand, more general spectra, which lead 
to different mass splitting, will be more or less sensitive to ISR. Figure 17 shows the effect of the 
mass splitting for qq events that pass the Jets+mHT event selection. Only the decay q — > qNi is 
allowed. The upper plot, with a fixed mass splitting of 100 GeV, shows a strong dependence on ISR. 
The lower plot, where is fixed at 100 GeV, shows a more balanced dependence on ISR and FSR. 
Some regions of Fig. 16 which overlap with the parameters in the lower part of Fig. 17 have g cascade 
decays, such as g — > qN2{-+ qqNi), and thus less dependence on in ISR. Cascade decays, however, 
are more likely to produce at least one lepton in the final state, leading to a loss of acceptance with a 
leptonic veto. 

The issue of the impact of ISR and FSR on high-p^ searches needs to be studied carefully. To date, 
some basic studies have been performed [22-24], but these do not address theoretical uncertainties. 
Currently, the systematic uncertainties assessed on the experimental limits, by increasing or decreasing 
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Jets+mHT: Fraction of Events with ISR 
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Figure 16. The fraction of events that pass the Jets+mHT selection with an additional jet from ISR in the 
cMSSM plane (mo, mi) for tan/3 = 10. 

the amount of radiation, are conservative. A better description could allow us to set more stringent 
limits on new particle production, or can be leveraged as a tool for discovery. Data driven techniques 
would be beneficial, perhaps based on observables from tt production (particularly pT{t + t)). 

6 Conclusions 

The textbook signature for SUSY at hadron colliders is Jets+fx, arising from the production and 
decay of the gluino g and the squarks q. An analysis of the limits arising from hadronic SUSY searches 
by CMS, and interpreted in the cMSSM, shows the current limits arc mainly excluding direct and qg 
production and decay. The exclusion of qq is an indirect exclusion of the the gluino 5 which mediates 
the process. Particle-level studies that follow the CMS event selections for the Jets+mHT and ut 
analyses, show that the dominant processes at the cross selection level are also the dominant ones after 
cuts. Our study also quantified the overlap between the two different selections (20-50%, depending 
on the process) and the impact of the lepton veto (up to a 50% loss of events). 

The importance of qq and qg is interesting, since interpretations of the data beyond the cMSSM 
often assume qq* and gg production. We suggest these new topologies and ranges of parameters for 
them. These processes need also to be considered for a construction of a cMSSM point in terms of 
simplified models. 
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Figure 17. The acceptance of events passing the Jets+mHT event selection divided by their dependence on 
an ISR jet, for a fixed mass and varying mass spHtting. 

We also recognize the large potential of electroweak gaugino production for exclusion or discovery of 
SUSY. Beyond the current limits, chargino pair and neutralino-chargino production are the dominant 
processes. These processes also impact the Icptonic searches for SUSY. In fact, the current CMS 
leptonic limits are bordering the region of the cMSSM plane (for the tan/3 = 10 slice) where the 
second heaviest neutralino decays to a Higgs boson and the lightest neutralino. 

Because of the relevance of qq production to the current limits, we studied in some detail the 
importance of initial state radiation for passing the SUSY search selections. For the cMSSM points 
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considered, final state radiation is more relevant, except for gluino masses less than 300 — 400 GeV. 



7 Appendix 
Squark-Gluino Reference Plot 

To ease the interpretation of figures displayed in the {mo, mi) plane, we collect here reference plots 
displaying contours of physical masses. 
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Figure 18. For reference, contours of squark and gluino masses in the cMSSM for tan/3 — 10, = 
GeV, sgn(ij.) > 0. A line with slope m — 0.9 has been superimposed to show the approximation correlation 
between Mq and Mg. 
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Figure 19. For reference, contours of sbottom, stop, and squark masses in the cMSSM for tan /3 = 10, Aq = 
GeV, sgn{fj,) > and tan ^ = 40, Ao = -500 GeV, sgn{^) > 0. 



Dominant Processes for tan/3 ~ 40 

We present a similar plot as in Fig. 2, but for a cMSSM point with a large value of tan/3. 
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Fractional Size:tanp=40Ao=-500 GeV,n > 
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Figure 20. For reference, the fractional size of various subprocesses as a function of {mo, mi) for tan/3 
40, Ao = -500 GeV, sgn{fj,) > 0. 



Squark-Squark and Squark-Antisquark Differences 

We provide additional information on the kinematic differences between and qq* production, namely 
the distribution of the rapidity y of the q{q*)- 

Pseudo-Analysis Details 

Here, we provide details on the tools used to simulate the detailed analyses performed by CMS, namely 
Jets+mHT and ax- For each value of (mg, mi) considered, the Pythia6[25, 26] Monte Carlo was used 
to simulate the production, decay, and parton showering of sparticles. The pT-ordered shower was 
used, with the underlying-event and hadronization effects turned off. Information about the progeny 
of jets was based on Monte Carlo truth information about the highest pt parton inside the jet. Jets 
were defined using the package Fast Jet [27, 28]. 

The cuts applied are based on the description in the CMS analysis studies. The lepton veto is 
applied, again, using Monte Carlo truth information and leptons from hard decays. 
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Figure 21. Comparison of tlie rapidity y of tlie produced q, q* for qq versus gg* production. 
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